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Recent fluorescence spectroscopy measurements of the turnover time distribution of single-enzyme 
turnover kinetics of /1-galactosidase provide evidence of Michaelis-Menten kinetics at low substrate 
concentration. However, at high substrate concentrations, the dimensionless variance of the turnover 
time distribution shows systematic deviations from the Michaelis-Menten prediction. This difference 
is attributed to conformational fluctuations in both the enzyme and the enzyme-substrate complex 
and to the possibility of both parallel and off-pathway kinetics. Here, we use the chemical master 
equation to model the kinetics of a single fluctuating enzyme that can yield a product through 
either parallel or off-pathway mechanisms. An exact expression is obtained for the turnover time 
distribution from which the mean turnover time and randomness parameters are calculated. The 
parallel and off-pathway mechanisms yield strikingly different dependences of the mean turnover 
time and the randomness parameter on the substrate concentration. In the parallel mechanism, the 
distinct contributions of enzyme and enzyme-substrate fluctuations are clearly discerned from the 
variation of the randomness parameter with substrate concentration. From these general results we 
conclude that an off-pathway mechanism, with substantial enzyme-substrate fluctuations, is needed 
to rationalize the experimental findings of single-enzyme turnover kinetics of /3-galactosidase. 


I. INTRODUCTION 
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Several biological processes that sustain life depend crucially on the catalytic activity of enzymes. A simple reaction 
mechanism for enzyme-catalysed reactions was first proposed in 1913 by Michaelis and Menten^ following the work of 
Wurtz, Henri and several others.^ According to the Michaelis-Menten (MM) mechanism, enzyme E binds with substrate 
S to form an enzyme-substrate complex ES^ which either dissociates irreversibly to form product P, regenerating the 
free enzyme E or dissociate reversibly to release the substrate: 


E + S%ES%E^ 

fc_i 


P 


( 1 ) 


In the above mechanism, ki and fc_i are the forward and backward rate constants for reversible substrate binding and 
substrate release steps respectively and ^2 is the rate constant for irreversible product formation and enzyme regeneration 
step. 

The classical MM kinetics, based on deterministic mass action kinetics, estimates the temporal variation of the concentra¬ 
tions of enzyme, enzyme-substrate complex and product once the initial concentrations are known. xhe deterministic 
nature of kinetics implicitly assume that all enzymes react at the same time to form products in bulk. In the steady-state, 
therefore, the initial rate of product formation is given by the classical Michaelis-Menten (MM) equation. 


2[E]o[S] 
[5] + Km ’ 


where v is the enzyme velocity. In the above expression, [Ajo = [E] + [ES] is the initial enzyme concentration and 
Km = {k-i + k 2 )/ki is the Michaelis constant. When the scaled reciprocal rate [E\q/v is plotted against the reciprocal 
concentration I/[S'], the MM equation yields a linear curve with the slope and intercept given by KM/k 2 and 1/^2 
respectively. This transformation, suggested by Lineweaver and Burk^ in 1934, is widely used to obtain the kinetic 
parameters of enzyme-catalyzed biochemical reactions. 

A crucial step in the Michaelis-Menten mechanism is the binding of the enzyme with a substrate to from an enzyme- 
substrate complex, dissociation of which leads to the formation of product. An enzymatic conformational state is, 
however, not fixed but undergoes incessant fluctuations due to random collisions with the surrounding solvent molecules. 
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The latter impart enough thermal energy to the enzyme to pass through the energetically unfavourable transition states 
that separate different conformational states. In an enzyme catalyzed reaction, therefore, different conformational states 
of an enzyme or an enzyme-substrate complex represent different kinetic intermediates which can interconvert between 
each other on time scales slower or faster than the product formation (or enzyme turnover) time.®“® How does the 
rate of enzymatic conformational fluctuations influence the turnover time of a single enzyme? Does the MM kinetics 
of a single fluctuating enzyme shows deviation from the classical MM equation? Do these fluctuations about the mean 
turnover time represent irrelevant random noise or carry useful information about the mechanistic pathway that can be 
characterized quantitatively? These are some of the questions that have just begun to be addressed in several recent 
studies on single-enzyme catalysis.®”^® 

In a recent single-enzyme turnover experiment based on fluorescence spectroscopy,®’® the turnover time for the repeated 
turnover of /3-galactosidase has been measured over a long time interval. From the frequency of occurrence of turnover 
times, the turnover time distribution of the enzymatic turnovers has been obtained, the first moment of which yields 
the mean turnover time. The turnover time distribution is mono-exponential at low substrate concentration and multi¬ 
exponential at high substrate concentration. Remarkably, the reciprocal of the mean turnover time exactly recovers the 
classical MM equation at low substrate concentration, 1/ (t) = implying that the mean turnover time obtained 

from the single-enzyme turnover time distribution is related to the (ensemble-average) enzymatic velocity in the classical 
MM kinetics. At high substrate concentration, however, the turnover time distribution shows multi-exponential decay 
resulting in systematic deviations from the MM kinetics. In particular, the dimensionless variance, which is expected 
to be unity for (single-pathway) MM kinetics is found to be greater than unity at high substrate concentration. This 
difference is attributed to conformational fluctuations in both the enzyme and the enzyme-substrate complex and to the 
possibility of both parallel and off-pathway MM kinetics. 

In the past few years, the randomness parameter has found widespread use as a new kinetic parameter to characterize 
temporal fluctuations in single-enzyme kinetics.® The randomness parameter, which is the dimensionless ratio of the 
variance and the square of the mean r = -——, provides a statistical measure of temporal fluctuations in the 
underlying mechanism in terms of the number of kinetic intermediates in the reaction pathway and their connectivity 
with respect to each other.Interestingly, a reaction mechanism with n sequentially connected kinetic states, I 

- A-n-^n-l-l, where occupancy time for each kinetic state is assumed to be exponentially distributed 
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uj{t) = ke~^*, always yields r = 1/n for going from 1 to n+1 state.Here, n represents the number of rate determining 
steps in a sequential mechanism with the same rate constant k. Thus, r = 1/n is the minimum amount of randomness 
that can be captured by the randomness parameter in a sequential mechanism.® This implies that the values of r < 1 
or r = 1 signify multiple or single rate determining step(s) respectively. The opposite limit of r > 1, on the other 
hand, implies off-pathway states as in inhibited enzyme mechanism.^® The randomness parameter can become greater 
than unity even when product formation occurs via parallel-pathway mechanism due to enzymatic interconversion. The 
latter mechanism has been proposed to rationalize the results of a recent single-enzyme turnover experiment on /3- 
galactosidase, where the randomness parameter increases from unity to greater than unity with the increase in substrate 
concentration.®’® At high substrate concentration, it saturates to a constant value greater than unity. 

The following two-state parallel-pathway MM mechanism serves as a minimal model to observe the effects of enzymatic 
conformational fluctuations in the MM kinetics:®’® 

Ex+S% ESx^ El + P, El ^ El 

fc-ii 

a] [a P] t/3 7l ll 

E2+S ES2 ^ El+ P, El ^ E2 (3) 

fc -12 

In this mechanism, a single enzyme, enzyme-substrate intermediate and the regenerated enzyme at any time t can exist 
in any one of the two conformer states given by Ei or E 2 ; ESi or ES 2 and A® or A® respectively. Here, a, (3 and 7 are 
the rate constants specifying interconversion between two enzyme, enzyme-substrate and regenerated enzyme conformers 
respectively. 

In a recent work,^^ an approximate expression for the turnover time distribution for the above mechanism has been 
obtained in terms of the individual turnover time distributions of the first Jei (t) and second fE 2 (t) enzyme conformers, 

fit) = WlfE^it) +W2fE2it), (4) 

where /ei (t) and second fE 2 (t) are the turnover time distributions for the independent (single-pathway) MM reaction 
steps, given by Ei ^ ESi El + P, El -5- Ei and E 2 ES 2 —>■ A® -I- P, A® —> E 2 respectively. Also, rci and W 2 
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are the steady-state probabilities for the enzyme to exist in either Ei or E 2 state. This form of expression includes 
the details of enzyme and enzyme-substrate conformational fluctuations in wi and W 2 . The above distribution shows 
mono-exponential decay at low substrate concentrations and multi-exponential decay at high substrate concentrations 
which is in qualitative agreement with the experimental findings. It does not, however, recover the Michaelis-Menten 
equation [Eq. Q] in the absence of enzymatic conformational fluctuations (a = /3 = 0). In the absence of this limiting 
behaviour, the conditions under which enzymatic conformational fluctuations influence the MM kinetics can not be fully 
characterized. 

In this work, starting from the chemical master equation for a single fluctuating enzyme following the parallel or 
off-pathway mechanisms, we obtain an exact expression for the turnover time distribution, the mean turnover time 
and the randomness parameter. The mean turnover time and randomness parameter exactly recover the MM kinetics 
in the absence of enzymatic fluctuations. In the presence of these fluctuations, however, the parallel and off-pathway 
mechanisms show strikingly different variations of the mean turnover time, and the randomness parameter, with substrate 
concentration. Further, the individual contributions of enzyme and enzyme-substrate fluctuations in the parallel-pathway 
mechanism are clearly discerned from the variation of the randomness parameter with substrate concentration. From 
our general analysis, we conclude that off-pathway mechanism with single product formation step in the sole presence 
of enzyme-substrate fluctuations is required to rationalize the results of the single-enzyme turnover experiment on 
/3-galactosidase. 

This paper is organized as follows. Section II describes the general features of the parallel-pathway MM mechanism 
with enzymatic conformational fluctuations. It also outlines the key steps of the calculation. Subsections A, B and C 
of Section III provide numerical analysis of the analytical expressions in terms of the turnover time distribution, the 
mean turnover time and randomness parameter for single-pathway, parallel-pathway and off-pathway MM mechanisms 
respectively. Section IV presents a brief summary of the results along with conclusions. The calculations are detailed 
in Appendices A and B. 

II. SINGLE-ENZYME MM KINETICS IN THE PRESENCE OF ENZYMATIC FLUCTUATIONS 


At the single-enzyme level, products are not formed in bulk but one at a time, resulting in a single enzyme turnover 
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in each catalytic cycle. The inherent stochasticity of a chemical reaction ensures that the repeated turnover of a single 


enzyme observed over a long time interval yields a distribution of turnover (or dwell) time.^^ For the parallel pathway 
MM mechanism [Eq. ([3])], the intrinsic stochasticity of an enzyme-catalyzed reaction and discrete integer jump in the 
number of reactants, intermediates and products at time t are accounted for in the chemical master equation approach 
(CME)^^“^® of stochastic processes.The effects of stochasticity are included by considering the number of enzymes 
enzyme-substrates {uESi), regenerated enzyme {n-Eo) and products (rip) as discrete random variables which can only 
take finite number of positive integral values as determined by stoichiometry. The CME for the parallel-pathway MM 
mechanism [Eq. (| 31 )] is given by 


dtP = 


— l) riEi + k-uC^El^EiS — ^)nEiS 

,i=l 

+ k 2 i(^EiS'^El'^E 0 — ^)nEiS^ P + a — l)n_Ei + P 

+/3 - l)n_EiS + (E^^^EejS - P 


( 5 ) 


where P represents the joint probability distribution given by P(nEi, riESi, ^^ 2 : nEo^np] t) and E is the step 

operator^® which operates on an arbitrary function f(x) yielding E/(a;) = f(x + l) and E“^/(a;) = f(x—l). This implies 
(E 3 ;E“^ — l)xf(x, y, z) = (x + i)f(x + l,y — 1, z) — xf(x, y, z). In the above equation k[^ = fcii[ 5 '] and k[2 = ^12 [S'] are 
the pseudo first-order rate constants. 

The CME is written for the time evolution of the joint probability of the number of each species involved in the chemical 
reaction.For a single enzyme, at any time f, different enzyme (or enzyme-substrate) conformational states are mu¬ 
tually exclusive implying PEi{t) = P(l, 0, 0,0, 0,0, 0; t), PESi(t) = 1^(0,1,0, 0,0, 0,0; t), PE 2 (t) = P(0, 0,1, 0,0, 0,0; t), 
Pss^ik) = ^*(0,0, 0,1, 0,0, 0; t).^® As a result, the chemical master equation reduces to the following set of coupled 




differential equations, given by 
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— + a)PEi (t) + k-iiPESi (t) + ckPe^ (t) 

KiPei (t) — (k-ii + k2i + P)PeSi (t) + PPes- 2 {t) 

— {k'i2 + a)PE2{P) + fc-i2^£;S2 W + oiPEi (t) 

^12^-E2 (t) + PPeSi {t) — {k-i2 + /3 + k22)PES2 W 

k2lPESi (0 + k22PES2 ■ ( 6 ) 

As the step ^ Ei is considered instantaneous, Pg O ~ 0 and does not appear in the above equations.Given 

this, the constraint, Psiit) + PESi{t) + Pe 2 {P) + PES 2 {t) = 1 has to be satisfied at all times. It is to be noted that the 
above differential equations can also be obtained directly by replacing the “concentration” descriptor in deterministic 
kinetics with the “probability” descriptor. However, in the context of enzymatic fluctuations, different conformational 
states only have probabilistic interpretation and a rigorous derivation of the above equations is possible only from a 
CME description of the kinetics. 

Since the enzyme turnover event is coupled to the product formation, the turnover time probability of an enzyme 
turnover to occur between time t and t + At is equal to the time for the product formation P in the same time interval. 
This implies f{t)At — APp{t), which in the limit of At —>• 0 yields the following expression for the turnover time 
distribution: 


dPEi (t) 
dt 

dPESi ji) 
dt 

dPE2 jt) 
dt 

dPES2 jt) 
dt 

dPp(t) 

dt 


f{t) = -= k2lPESi{t) + k22PES2{t) 


( 7 ) 


To obtain turnover time distribution from the above expression, Eqs. (ED can be solved to obtain PpSiit) and PES 2 {k)- 
In the next section, we provide a detailed analysis of the turnover time distribution and its statistical moments to 
understand the role of enzymatic conformational fluctuations in single-enzyme MM kinetics. 








III. TURNOVER STATISTICS OF A SINGLE ENZYME 


The turnover time or dwell time distribution and its statistical moments such as the mean (t) and the variance 
carry useful information about the mechanistic pathway of an enzyme catalysed reaction.®’^® In subsections (B) and (C) 
of this section, we analyse the statistical properties of an enzyme catalyzed reaction which follows the parallel-pathway 
and off-pathway mechanisms respectively. In subsection (A) below, we first recapitulate the results of a single-enzyme 
catalysed reaction following the single-pathway MM mechanism [Eq. (IT|)]. 

A. Turnover statistics of a single enzyme following the single-pathway MM mechanism 

The turnover time distribution of a single-enzyme catalysed reaction has been obtained previously from the CME 
approach,which yields 


fit) 


k2h[S] 

2A 


^(A-B)t _ ^-(A+B)t 


( 8 ) 


where A = sj (/ci[5'] -I- /c_i -I- ^ 2 )^ — 4A:iA;2[5']/2 and B = (fci[5'] -I- k-i k2)l2. From the above equation, the first and 
second moments are given by (t) = dt tf{t) and dt fit). These expressions yield the mean turnover 

time and the randomness parameter as^^’^^ 


it) 


-I- k-i + k2 

kik2[S] 


(9) 


and 


(f^) - (t)" _ (fcl[g]+fc-l)"+2fc2fc-l+fc| 
iki[S] + k-i + k2)'^ 


( 10 ) 


respectively. Comparison of Eq. 0 with Eq. ([5]) yields 1/ (t) = v/[E\q implying that the reciprocal of mean turnover 
time obtained from the single-enzyme turnover time distribution is related to the (ensemble-average) enzymatic velocity 
in the classical MM kinetics. The minimum value of r can be obtained using dr/d[S] |[s]=[s]' = 0 resulting in = 
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(^2 + /c_i)/fci. At this substrate concentration, = ^’(['5']') yields 

= 2(fc2 + fc-i)- 

In the limiting case of fc_i = 0, the minimum value of r is given by Tmin = 'f([>S']0 = 1/2 with = ^ 2 - This implies 

that there are two rate determining steps (n = Xjrrain = 2) corresponding to the substrate binding and the product 
formation steps respectively. The above equations will be used to analyse the limiting behaviour of the MM kinetics in 
the presence of enzymatic fluctuations. 

B. Turnover statistics of a single fluctuating enzyme following the parallel-pathway MM mechanism 

1. Turnover time distribution 

The turnover time distribution for the parallel-pathway MM [Eq. Q] mechanism is given by Eq. 0. To obtain an exact 
expression for the turnover time distribution, we first solve the coupled differential equations [Eq. ([6])] in the Laplace 
domain, the details of which are presented in Appendix A. Since the expression for the turnover time distribution in the 
time domain is unwieldy [Eqs. (IA4I) and (IA5I) ]. we directly evaluate the inverse Laplace transform numerically to obtain 

fit)- 

Fig. (1) shows the temporal variation of the turnover time distribution for the parameter values reported in the figure 
caption. The turnover statistics of a single fluctuating enzyme is mainly governed by two distinct time scales - the 
time scale of interconversion between enzymatic conformational states and the time scale of product formation which is 
linked to the enzyme turnover event. In what follows, we refer to enzymatic fluctuations as the ones that originate from 
interconversion between enzyme [Ei ^ A 2 ] and enzyme-substrate [ESi ^ ES 2 \ conformational states. The enzyme or 
enzyme-substrate conformational fluctuations respectively refer to the presence of either the former or the latter. 

Fig. (la) shows the effect of enzymatic conformational fluctuations on turnover time distribution for parallel-pathway 
mechanism. In the absence of interconversion between enzymatic conformers, a = /3 = 0, the decay is mono-exponential 
[brown curve]. The filled circles represent turnover time distribution [Eq. ([5])] for the MM mechanism [Eq. ([IJ] which 
exactly coincides with the brown curve indicating that single-pathway MM mechanism is followed in the absence of 
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FIG. 1: Single-enzyme turnover time distribution as a function of turnover time for the parallel-pathway MM mechanism (in 
non-dimensional units) for (a) different rates of enzyme and enzyme-substrate interconversion characterized by rate constants a 
and respectively and (b) asymmetry in the rate constants for substrate binding, substrate release and product formation steps. 
Common parameter values in (a) and (b) are fcii = fc_ii = ~ ^-12 ~ 0.5, fc 2 i = 5, ~ 0.1, a = j3 = 1*-“^ and [S'] = 50. 

The parameter values which were varied in subfigures (a) and (b) are indicated as superscripts and the corresponding values are 
indicated as subfigure legends. Filled circles correspond to the turnover time distribution [Eq. (|8])] for the (single-pathway) MM 
mechanism [Eq. ([T])] with fci = fe_i = 0.5, fc 2 = 5 at [S] = 50. The latter exactly coincides with the turnover time distribution for 
the parallel-pathway MM mechanism (a) in the absence of enzymatic fluctuations a = /3 = 0 (brown curve), and (b) symmetry 
in the catalytic rate constants ^21 = ^22 (red curve). 


enzymatic fluctuations. The mono-exponential decay is also observed for the case when enzymatic fluctuations occur 
faster on the time scale of the catalytic step [green curve]. Single time scale of decay implies that the effects of fast 
enzymatic fluctuations are averaged out on the time scale of product formation. 

In between these limits, the decay profiles are multi-exponential implying multiple competing reaction pathways. For 
the symmetrical case, a = j3, multiexponential decay is observed when enzymatic conformational fluctuations occur 
slower on the time scale of the catalytic step. For the asymmetrical case, a ^ P, multiexponential decay is observed 
when the time scale of either enzyme or enzyme-substrate interconversion is slower than the catalytic step. All these 
curves satisfy the condition, ^21 ^ ^ 22 - Fig- (lb) exemplifies the latter condition for the symmetrical case, a = P = 1. 
When enzymatic conformational fluctuations occur slower on the time scale of the the catalytic step, then asymmetry 
in the rate constants for either the catalytic steps A; 2 i yf /C 22 or the substrate binding steps kn ^ ki 2 is necessary to 
observe multiexponential decay. On the contrary, asymmetry in the substrate release step fc_ii ^ k -12 shows mono¬ 
exponential decay [black curve]. The same is also true in the absence of asymmetry in the catalytic steps, which results 
in mono-exponential decay [red curve] implying that single-pathway MM mechanism is followed [filled circles]. 
Although the present description only accounts for fluctuations between two (instead of n) enzyme or enzyme-substrate 
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FIG. 2: Temporal variation of the turnover time distribution for different substrate concentration for the parallel-pathway MM 
mechanism (in non-dimensional units). The turnover time distribution shows monexponential decay at low substrate concentration 
and multiexponential decay at high concentration. The parameter values are fcn = fc_ii = fci 2 = k -12 — 0.5, ^21 = 5, k 22 = 1, 
and a = P = 1. Filled circles correspond to the turnover time distribution [Eq. ([SI] for the (single-pathway) MM mechanism 
[Eq. O] with fci = k-i = 0.5, fc 2 = 5 at [5] = 1. The latter shows deviation from the turnover time distribution for the 
parallel-pathway MM mechanism at the same substrate concentration (blue curve). 


conformers, it is sufficient to capture the conditions under which the presence of enzymatic fluctuations can result in 
multi-exponential decay. Fig. (2), for instance, captures the salient features of the turnover time distribution as a 
function of substrate concentration. This include mono-exponential decay at low substrate concentration and multi¬ 
exponential decay at high substrate concentration, in qualitative agreement with the single-enzyme turnover experiment 
on /3-galactosidase. At low substrate concentration, however, the turnover time distribution for the parallel-pathway 
MM mechanism [Eq. ([7])] does not coincide with the turnover time distribution for the (single-pathway) MM mechanism 
[Eq. (|8])]. This indicates that the MM behaviour is not recovered at low substrate concentration in spite of the mono¬ 
exponential nature of the decay curve. The consequences of the latter are discussed below. 

It is to be noted that in the context of two-state model of enzymatic conformational fluctuations presented here, the 
asymmetry in the rate constants for the catalytic (or substrate binding) step amounts to dynamic disorder in the n-state 
model implying a broad distribution of ^2 (or fci).^^ 


2. First statistical moment and mean turnover time 


We can obtain more statistical information from f{t) by taking its first moment which yields the mean turnover time, 
(t) = tf(t)dt. In the Laplace domain, it can be obtained from (t) = — , where F{s) is the Laplace transform 


s=0 










of f{t), given by F(s) = k 2 iPESi{s) + fc 22 ^’£;S 2 («)■ This yields 


12 


it) 


As — k2iAi — A:22^2 

A 4 

£;[5]2 + F[5] + G 
id[^]2 + I[S] 


( 12 ) 


where, 

E = kiiki2{k22 + 2f3) 

F = A:ii(afc_i2 + Q;fc22 + 2a/3 + A:_i2/3) 

+ ^ 12 (^- 11^22 + Pk-ii + ^ 21^22 + k 2 lP 
+ ^22/3 + Q;fc_ii + Oik 21 + 2ci/3) 

G = ‘2oi{k—\xk—\2 + k—iik22 + Pk—w + k—\2k2i 

-\-k21k22 + Pk2i + 13 k—12 + / 3 k 22 ) 

H = kiiki2{k2ik22 + k2i(3 + k22f3) 

I = aki2{k-iik22 + ^21^22 + Pk2l + /3fc22) 

+ Q;fcii(fc2lfc-12 + felfc22 + /3fc21 + /3fe2)- 


The enzymatic velocity of single fluctuating enzyme can be obtained from the reciprocal of the mean turnover time, 


1 H[Sf + I[S] 

(t) ~ £;[5]2+F[5] + G’ 


(13) 


From the above expression it can be easily seen that the Michelis-Menten equation, which is characteristically hyperbolic 
in the substrate concentration [S], does not hold in the presence of enzymatic conformational fluctuations. However, 
there are certain limiting conditions under which the above equation recovers the Michaelis-Menten equation or the 
Michaelis-Menten like equation : 

(i) Of = /3 = 0 exactly recovers the Michaelis-Menten equation, v = , 
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where H' = /ciifc 2 i, E' = kn and F' = k-n + A: 2 i- 

(ii) a = 0 yields the Michaelis-Menten like equation given by n = , 

where F' = ^[fciifc_i 2 + A:i 2 (/c-ii + ^21 + ^ 22 )] + ki 2 {k-iik 22 + ^ 21 ^ 2 ) 

(iii) a = 0 and /3 —)• cxd yields the Michaelis-Menten like equation given by u = e^s\+F' ^ 
where E' = 2fciifci2, F' = fciifc_i 2 -I- fci 2 (fc_ii -I- fei -I- ^ 22 ) and H' = fcii/ci 2 (fc 2 i -I- k 22 ) 

(iv) a ^ 00 yields the Michaelis-Menten like equation given by u = , 

where F' = fcii(fc_i 2 -I- ^2 + 2/3) -I- fci 2 (/c_ii -I- k 2 i + 2/3), G' = G/a and /' = I/a. 

(v) a 00 and /3 —)■ 00 yields the Michaelis-Menten like equation given by w = g/jgj^g, , 
where F' = 2(fcii -|- ^ 12 ), G' = 2(A:_ii -|- k -12 + k 2 i + ^ 22 ) and /' = (^21 -I- /c 22 )(fci 2 + fcii)- 

The above analysis shows that the MM equation is exactly recovered only in the absence of enzymatic conformational 
fluctuations [condition (i)]. In order to analyse Eq. (jl2p in greater detail, we numerically evaluate the mean turnover time 
(t) as a function of the reciprocal of the substrate concentration 1/[<S'], which yields the single-enzyme Lineweaver-Burk 
(LB) curves [Fig. (3)]. 

Fig. (3a) graphically depicts the conditions under which the MM or MM like behaviour can be obtained. In the absence 
of enzymatic conformational fluctuations, a = /3 = 0 [condition (i)], the MM equation is exactly recovered yielding a 
linear curve. This is shown as the collapse of the linear red curve with filled circles, obtained from the variation of 
[E]q/v with respect to l/)^] in the MM equation [Eq. ([5])]. The conditions (ii)-(v), represented by green, black, blue 
and pink curves respectively yield linear curves with intercepts and slopes different from that of the MM equation [red 
curve], thereby showing the MM-like behaviour. For the symmetrical case, a = /3, the non-linearity in the LB curve 
occurs when enzymatic fluctuations occur on time scale longer than the catalytic step [brown curve]. All these curves 
satisfy the condition, k 2 i ^ ^ 22 - When ^21 = ^ 22 ) all linear and non-linear curves merge with the red curve to exactly 
recover the MM equation implying that single-pathway MM mechanism is followed under these conditions [Fig. (3b)]. 
Fig. (3c) shows the effect of enzyme or enzyme substrate conformational fluctuations on the mean turnover time [con¬ 
ditions (ii)-(iv)]. In the limiting cases of a = 0, ^ 00 and a —>■ 00 , the MM-like behaviour is recovered corresponding 

to the conditions (iii) and (iv) respectively. The sole presence of enzyme-substrate conformational fluctuations [merged 
brown and purple curves] also yields MM-like equation [condition (ii)]. The sole presence of enzyme conformational 
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1/[S] 



1/rsi 


FIG. 3: Single-enzyme Lineweaver-Burk plots for (t) versus 1/[S'] for the parallel pathway MM mechanism (in non-dimensional 
units) for (a) different rates of enzymatic interconversion characterized by the rate constants a and for fe 2 i yf fc 22 , (b) same as 
(a) but in the absence of asymmetry in the rate constants for the catalytic steps fc 2 i = fc 22 , (c) different rates of conformational 
fluctuations arising solely from enzyme or enzyme-substrate interconversion characterized by the rate constants a ox P respectively 
and (d) asymmetry in the rate constants for substrate binding, substrate release or catalytic steps. Common parameter values 
in (a)-(d) are fen = fc_ii = fe^^ = ^^-12 ~ 0-5, fe 2 i = 5, ~ 0.1 and a = P = 'pjjg parameter values which were 

varied in subhgures (a)-(d) are indicated as superscripts and the corresponding values are indicated as subfigure legends. For the 
chosen parameter values in (a), the green, black and pink curves almost merge onto each other. Filled circles correspond to the 
mean turnover time [Eq. ([91] for the (single-pathway) MM mechanism [Eq. ([T])] with fei = fe_i = 0.5, fe 2 = 5. 


fluctuations, on the other hand, yields non-linearity in the LB curves even when enzyme conformational fluctuations 
occur faster on the time scale of the product formation step [blue, pink and gray curves]. Irrespective of the time scale of 
enzyme conformational fluctuations [blue, pink and gray curves], the parallel-pathway MM mechanism approaches the 
single-pathway MM mechanism [red curve] at high substrate concentration. In the sole presence of enzyme-substrate 
conformational fluctuations, however, the brown and purple curves do not approach the red curve implying that the 
parallel-pathway MM mechanism is followed at all substrate concentration. Once again, all these curves satisfy the 
condition fe 2 i 7 ^ ^ 22 - 

To explore the condition of asymmetry in the rate constants further, Fig. (3d) shows the LB curves for the symmetrical 
case, a = P = 1, when enzymatic conformational fluctuations occur slower on the time scale of the catalytic step. 
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For this condition, asymmetry in the rate constants for the catalytic step shows nonlinearity in the LB curve [purple 


curve]. Asymmetry in the substrate binding step fcn ^ fci2 shows linear MM- like behaviour with a shift from parallel 
to single-pathway mechanism at high substrate concentration [green curve]. Asymmetry in the substrate release step 
k_ii 7^ fc_i2, however, shows close agreement with the MM equation [blue curve] implying that single-pathway MM 
mechanism is followed. 

Figs. (3a), (3c) and (3d) show that the mean turnover time captures the substrate concentration dependent switch 
from parallel to single-pathway MM mechanism in terms of the approach of the non-linear LB curves or linear MM-like 
curves towards the exact MM equation [Eq. ([2])] at high substrate concentration. In next section, we show how the 
randomness parameter corroborate these results. 


3. Second statistical moment and randomness parameter 


The randomness parameter (or squared coefficient of variation)^®’^'^ is a dimensionless ratio of the variance to the mean- 
square, r = quantifies the magnitude of temporal fluctuations in a given reaction pathway which can be 

correlated with the underlying reaction mechanism. The values of r = 1 or r < 1 suggest a sequential mechanism with 
one or more than one rate determining steps respectively. The opposite limit of r > 1 can arise if there are off-pathway 
states with single pathway for product formation or parallel-pathway states with more than one pathway for product 
formation. 

To obtain r = , we first evaluate the second moment of /(t), which is given by (t^) = f{t)dt or (t^) = 

, where F(s) = k 2 iPESiis) -I- k 22 PES 2 {s)- This yields 


d^F{s) 

ds^ 


s=0 


2\ ^11^21 ~ 2A2 2A3(A:21^1 + ^ 22 ^ 2 ) 2A3 


r) = 


A4 


A2 


A? ’ 


(14) 


where the expressions for A 2 , A 3 , A 4 , Ai and A 2 are given in Appendix A. Eqs. (fT^ and (fTT)) are then used to obtain 
r. The variation of the randomness parameter with substrate concentration is shown in Figs. (4a)-(4d). A common 
feature of these figures is that r = 1 at low substrate concentration as substrate binding is the rate determining step. 
Fig. (4a) shows the effect of asymmetry in the rate constants for the catalytic step when enzymatic conformational 
fluctuations occur slower on the time scale of the catalytic step. For perfectly symmetrical case, k 2 i = k 22 , the substrate 









FIG. 4: Randomness parameter as a function of substrate concentration (in non-dimensional units) for parallel pathway MM 
mechanism for (a) asymmetry in the rate constants for product formation step; (b) for asymmetry in the rate constants for substrate 
binding and substrate release steps; (c) different enzymatic and enzyme-substrate interconversion rate constants characterized 
by a = /? and a = 0 , /I 7 ^ 0 respectively and (d) different enzymatic and enzyme interconversion rate constants characterized 
a ^ /3 and a 7 ^ 0,/3 = 0 respectively. Common parameter values in (a)-(d) are fcn = fc_ii = = 0.5, k 2 i = 5 and 

^22 ~ 0.1 and a = P = fChC) parameter values which were varied in subfigures (a)-(d) are indicated as superscripts and the 

corresponding values are indicated as subfigure legends. Dashed black curves in (a)-(c) correspond to the randomness parameter 
[Eq. (1101) ] for the (single-pathway) MM mechanism [Eq. ([T}] with fci = fe_i = 0.5, k 2 = 5. 


concentration dependence of the randomness parameter shows initial decrease from unity followed by convergence to 
unity at high substrate concentration. The latter suggests that the single-pathway MM mechanism [dashed black curve 
obtained from Eq. (ITUl) ] is followed at all substrate concentration [purple curve]. Following the analysis of Eq. (fTT|) . 
it can be readily seen that the randomness parameter attains the minimum value of = 0.545 at [S']' = 11 for the 
parameter values fc_ii = fc_i 2 = /c_i = 0.5 and ^21 = ^22 = /c 2 = 5. At high substrate concentration, the randomness 
parameter approaches unity as the product formation step is the rate determining step. For the asymmetrical case, 
^21 ^ ^ 22 ) the randomness parameter first increases from its initial value of unity before saturating to a constant value 


greater than unity at high substrate concentration. More the asymmetry in the catalytic rate constants, higher the 
value r at which the saturation occurs. This suggests substrate concentration dependent switch from single-pathway to 
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parallel-pathway MM mechanism for the condition k 2 i ^ ^ 22 - 

The effect of asymmetry in the rate constants for the substrate binding and the substrate release steps are depicted 
in Fig. (4b) for the case where enzymatic fluctuations occur slower on the time scale of the catalytic step. Unlike the 
case of asymmetry in the catalytic rate constants discussed above, asymmetry in the substrate binding rate constants, 
fell 7 ^ ki 2 , shows that the randomness parameter instead of saturating to a constant value greater than unity, decreases 
to unity at high substrate concentration. The latter implies that the parallel-pathway MM mechanism switches to single¬ 
pathway MM mechanism at high substrate concentration. Asymmetry in the substrate release step k-n 7 ^ k -12 [purple 
curve], almost coincides with the single-pathway MM mechanism [dashed curve] implying that the single-pathway MM 
mechanism is followed at all substrate concentration. 

Fig. (4c) shows the effect of enzymatic conformational fluctuations for the symmetrical case, a = j3. Single-pathway 
MM mechanism is followed in the absence of enzymatic fluctuations [a = /3 = Oj. The same is also true for the case 
when the enzymatic conformational fluctuations occur faster on the time scale of the catalytic step [blue curve]. As the 
time scale of enzymatic conformational fluctuations become slower than the catalytic step, the randomness parameter 
first increases from its initial value of unity and then saturates to a constant value greater than unity at high substrate 
concentration. Slower the time scale of enzymatic fluctuations compared to the catalytic step, more the randomness. 
These results imply a shift from single-pathway MM mechanism at low substrate concentration to parallel-pathway MM 
mechanism at high substrate concentration, whenever enzymatic fluctuations [a = /3] occur slower on the time scale of 
the catalytic step. Interestingly, the sole presence of enzyme-substrate conformational fluctuations (a = 0) shows similar 
trends [dashed curves]. The latter suggest that irrespective of the time scale of enzyme conformational fluctuations, 
the switch from single to parallel-pathway MM mechanism at high substrate concentration is mainly determined by the 
competition between the time scales of enzyme-substrate interconversion and product formation steps. 

The asymmetrical case oi a ^ j3 shows that the distinct contributions of enzyme (/3 = 0) and enzyme-substrate {a = 0) 
conformational fluctuations can be clearly discerned from the variation of the randomness parameter with substrate 
concentration. Fig. (4d) shows that the sole presence of enzyme conformational fluctuations (/3 = 0) results in giant 
temporal fluctuations, characterized by very high noise-to-signal ratio (r ^ 1 ), even when enzyme conformational 
fluctuations occur faster on the time scale of the catalytic step. At very high substrate concentration, the randomness 
parameter becomes unity as the product formation becomes the rate determining step, implying a shift from parallel to 
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FIG. 5: Single-enzyme turnover statistics for the off-pathway MM mechanism [Eq. (IlSIl l (in non-dimensional units), (a) Turnover 
time distribution as a function of time at different substrate concentrations, (b) Single-enzyme Lineweaver-Burk plot based on 
the off-pathway MM mechanism [Eq. (1161) 1 shows that the MM equation is recovered at low substrate concentration, (c) Single¬ 
enzyme Lineweaver-Burk plot based on the off-pathway mechanism [Eq. dlBIl ] shows deviation from the MM equation at high 
substrate concentration, (d) The randomness parameter increases from r = 1 at low substrate concentration to r > 1 at high 
substrate concentration when enzyme-substrate interconversion rate is slower than the catalytic step. The parameter values for 
(a)-(d) are k\ = fe_i = 0.5, fe = 5 with /3 = 1 in (a). Filled circles in (a)-(c) and black curve in (d) correspond to the turnover 
time distribution, mean turnover time and randomness parameter for the (single-pathway) MM mechanism obtained from Eqs. 
m, ® and CQJ respectively with parameter values fei = fc_i = 0.5 and k 2 = 5. 


single-pathway MM mechanism. Interestingly, even a small increase in the value (3 [dashed curves in Fig. (4d)], suppresses 
the magnitude of these giant fluctuations to almost half its original value. Also, instead of decaying to unity at very 
high substrate concentration, the randomness parameter saturates to a constant value greater than unity, implying 
parallel-pathway mechanism is followed at high [S']. Thus, in the presence of enzymatic conformational fluctuations with 
a ^ (3, the randomness parameter exceeds unity when enzyme-substrate conformational fluctuation occur slower on the 
time scale of the catalytic step [dashed curves in Fig. (4d)]. In the opposite limit, r < 1 [dotted curves in Fig. (4d)]. 
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C. Turnover statistics of a single enzyme following the off-pathway MM mechanism 

/3-galactosidase is an enzyme that catalyzes the hydrolysis of lactose and is known to follow the MM kinetics. In a recent 
experiment based on fluorescence spectroscopy, the turnover time for the repeated turnover of a single /3-galactosidase 
was observed over a long time interval from which the turnover time distribution was obtained.® The latter showed 
mono-exponential decay at low substrate concentration and multi-exponential decay at high substrate concentration. 
The single-enzyme Lineweaver Burk plot, obtained from the variation of the mean turnover time with the reciprocal of 
the substrate concentration, showed exact agreement with the MM equation [Eq. ([2])] at low substrate concentration. 
The variation of the randomness parameter with substrate concentration, on the other hand, was found to be unity at 
low substrate concentration and greater than unity at high substrate concentration. 

The parallel-pathway MM mechanism, discussed in the previous section, can capture certain features of this experiment. 
This include (a) mono-exponential decay of the turnover time distribution at low substrate concentration and multi¬ 
exponential decay at high substrate concentration [Fig. (2)], and (b) the increase in the randomness parameter from 
unity to greater than unity with the increase in the substrate concentration [Figs. (4a) and (4c)]. Both these results 
require that enzymatic [or enzyme-substrate] conformational fluctuations are slower on the time scale of the catalytic 
step and k 2 i ^ k 22 - However, the same conditions that recover (a) and (b) can not recover the MM equation [Eq. 
@] at low substrate concentration. This effect is clearly captured in Fig. (2), which shows that in spite of the mono¬ 
exponential decay of the turnover time distribution at low substrate concentration [blue curve], the MM kinetics [filled 
circles corresponding to Eq. ®] is not recovered. 

Below, we show that experimental findings can be rationalized on the basis of the following off-pathway MM mechanism: 

Ei + S ^ ESi E° + P, E° % El 

k-i 

/ 31 113 

ES2, (15) 


where product formation occurs via single pathway mechanism. 

The turnover time distribution for the off-pathway mechanism has been evaluated in Appendix B. Fig. (5a) shows 
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the dependence of the turnover time distribution on substrate concentration. For enzyme-substrate conformational 
fluctuations occurring slower on the time scale of the catalytic step, the turnover time distribution is mono-exponential at 
low substrate concentration and multi-exponential at high substrate concentration. However, in contrast to the parallel¬ 
pathway mechanism, the turnover time distribution for the off-pathway mechanism at low substrate concentration [blue 
curve] exactly recovers the MM kinetics [filled circles corresponding to Eq. ([8])]. 

The first moment of the distribution yields 


_ 1 _ fcifc2[>5'] , . 

"" (t) 2fci[5] + (fc_i + fc2)’ ^ ^ 

which recovers the MM equation at low substrate concentration [Fig. (5b)], but shows deviation from it at high 
substrate concentration [Fig. (5c)]. With the increase in the substrate concentration, the randomness parameter shows 
an increase from its initial value of unity to greater than unity whenever enzyme-substrate conformational fluctuations 
are slower on the time scale of the catalytic step [Fig. (5d)]. The present analysis, thus, shows that an off-pathway MM 
mechanism, with substantial enzyme-substrate fluctuations, is required to rationalize the results of the single-enzyme 
turnover experiment on /3-galactosidase. 


IV. SUMMARY AND CONCLUSIONS 

In this work, we use the chemical master equation to model the kinetics of a single fluctuating enzyme that can yield 
a product through either parallel or off-pathway mechanisms. We obtain an exact expression for the turnover time 
distribution, from which the mean turnover time and randomness parameters are evaluated. 

For parallel-pathway MM mechanism, the reciprocal of the mean turnover time yields the exact MM equation in the ab¬ 
sence of enzymatic conformational fluctuations. The presence of enzyme, enzyme-substrate or enzymatic conformational 
fluctuations, in contrast, yield deviations from the MM equation whenever the rate constants for product formation (or 
substrate binding) steps are unequal. The variation of the randomness parameter with substrate concentration clearly 
discerns the distinct contributions of enzyme and enzyme-substrate fluctuations in the parallel-pathway mechanism. 
The randomness parameter is found to be greater than unity at high substrate concentration when enzyme-substrate 
conformational fluctuations are slow on the time scale of the catalytic step, implying that parallel pathway MM kinetics 
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is followed at this concentration. When only enzyme conformational fluctuations are allowed, in contrast, the random¬ 
ness parameter is greater than unity as long as enzyme conformational fluctuations occur slower on the time scale of 
the substrate binding step. At high substrate concentration, when the catalytic step becomes the rate determining step, 
the randomness parameter decays to unity implying the switch from parallel to single-pathway MM kinetics with the 
increase in the substrate concentration. 

The mean turnover time for the off-pathway MM mechanism, allowing only enzyme-substrate conformational fluctu¬ 
ations, yields the exact MM equation at low substrate concentration, but shows deviation from it at high substrate 
concentration. At these high concentrations, the randomness parameter is greater than unity whenever the enzyme- 
substrate conformational fluctuations are slow on the time scale of the catalytic step. 

From these general results, we conclude that an off-pathway mechanism with substantial enzyme-substrate conforma¬ 
tional fluctuations is needed to rationalize the results of the single-enzyme turnover experiment on /3-galactosidase. It 
is conceivable that the turnover statistics of different enzymes reveal different aspects of enzyme and enzyme-substrate 
conformational fluctuations in single-enzyme kinetics. The present analysis shows that the distinct contributions of en¬ 
zyme and enzyme-substrate conformational fluctuations in parallel or off-pathway mechanisms can be clearly discerned 
from the variation of the mean turnover time and the randomness parameter with substrate concentration. 



Appendix A: Exact calculation of the turnover time distribution for parallel-pathway MM mechanism 
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Eqs. (|ni) can be solved by taking the Laplace transform. Taking the Laplace transform and applying the initial conditions, 
-Pei(O) = 1, PeSi{0) = 0, PE 2 W = 0i PeS2{0) = 0 and Pp{0) = 0, the coupled differential equations reduce to the 
following set of algebraic equations: 

(s + k[^ + a)PEi (s) - fc-iiPeSi (s) - (s) = 1 

(s + fc_ii + P + k2i)PESi{s) — k[^PEi{s) — I 3 PeS 2 = 0 
(s + + a)PE2 (s) - k-i2PES2 (s) - olPe^ (s) = 0 

(s + /3 + fc_i2 + k 22 )PES 2 (s) - (s) “ PPpSi (s) =0 (Al) 

Rewriting the above equations in the matrix form yields 


(5 + + Of) 

—fc_ll 

— 0 

0 

\ 


Pei (s) 


fd 

— h' 

ft-ll 

(s + fc_ii + /3 + A: 2 i) 

0 

-p 



PeSi (s) 


0 

— 0 

0 

(s + k^i2 

—k-i2 



Pe2 (s) 


0 

V 0 

-/? 

— k' 

^12 

(s + /3 + fc_i2 

+ ^ 22 )j 


\PeS2 (s) y 


w 


Pei (s) 


(5 + k^w 0 ) 

-k-11 

—a 

„ 1 

-1 

fA 

PeSi (s) 


— k' 

^11 

(s + k-ii + P + k 2 i) 

0 



0 

Pe2 (s) 


—0 

0 

(s + ki 2 + a) 

—k-12 


0 

^PeS2 (s) y 


1 0 

-p 

— k' 

'>'12 

(s + /3 + k-i 2 + A 22 ) y 


[v 
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-Pei (s) 
PeSi (s) 
Pe2 (s) 

y-PeSa (s) j 


(g^-|-AQS^-|-.Bog-|-OQ) 
(s‘^+AiS^+A2s2+A3S + A4) 

-\-Ais-\-Bi) 

(s'^+AiS^+A2S^+A3S + A4) 

_ (q5^+.Do5+-E’o) _ 

(s^+Ais3+A2S^+A3S+A4) 

_ (A 25 +B 2 ) _ 

\ (s‘^+AiS^+A2S^+A3S+A4) j 


where, 


aIq — ^—11 H“ k—\2 ^21 H“ ^22 H“ 2/? H- /ci2[5'] H- g; 

^0 = (fc-11 + fc2l)(fc-12 + ^2) + (fc-11 + ^21 + k-i2 + ^22)/? 

+ (A:i 2[5'] + a)(fc_i2 + ^22 + 2/3 + k-n + ^ 21 ) ~ ^12^-12[S'] 

Co = (fci2[>S'] + a)[(k-ii + fc2i)(fc_i2 + A:22) + {k-ii + ^21 + A:_i2 + fc22)/3] 

-A:i2A:_i2(fc-ii + fc2i + /3)[S'] 

Do = ai^k—ii + k—i2 + ^21 + fe2 + 2/3) 

-P '0 = (^{k-ii + k 2 l){k-i 2 + k 22 ) + CK/3(fc-12 + P 22 + P -11 + P 21 ) + PiiP-12/5[>S'] 
Ai = Piifci 2 [S']^ + Pii(fc —12 + P 22 + Of + /3)[S'] 

Bi = Piifci2[5']^(p22 +/5) + Q;fcii[S'](fc_i 2 + P 22 +/3) + a/3fci2[S'] 

A 2 = {kiiP + ki2a)[S] 

B 2 = Piifci2/3[S']^ + fci2[5']a(fc_ii + ^21 +/3) + Piia/3[S'] 



and 
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Ai — {kii + A:i 2 )[<S'] + 2[a + /3) + k—n + fc_i 2 + fei + ^22 
A 2 = fcllfcl2[5']^ + fcli[S'](fc_i2 + ^21 + fcfc22 + a + 2/3) 

+ fcl 2 [s](fc —11 + ^21 + ^22 + Q; + 2/3) 

+ (fc_iiA:_i 2 + k-iik22 + fc-ii/3 + fc—12^21 + ^ 1^22 

+/c2i/3 + ^—12/3 + ^22/3 H—|- 2 Q;fc_ii + 2 q;/c 21 + 2oik —12 
+2ak22 + 4q;/3 

A 3 = fciifci 2 [S']^(A :21 + ^22 + 2/3) + fcii(fc 21^-12 + ^ 21^2 +/3fc21 
+A:_i 2/3 + A: 22/3 + q;/c 21 + otk —12 + cxk 22 + 2q;/3)[S'] 

+^ 12 (^- 11^22 + k-iiP + k2ik22 + ^ 21/3 + A:22/3 

+Q;fc_i2 + Q;fc2i + aA:22 + 2a/3)[S'] 

+ 2 (Q;fc_iifc_i 2 + ak-iik22 + ajSk-ii + ak-i2k2i 

-\-ak 21 k 22 + a(5k2i + a(3k-i2 + a(5k22) 

A 4 = ^ 11 ^ 12 (^ 21^2 + ^ 21/3 + A:22/3)[>S']^ 

+ Q;(fc_llfci 2 fc 22 + ^ 12 ^ 21^22 + ^12^21/3 + ^12^22/3 
+ ^11^21^-12 + ^11^21^22 + Pkiik2l + /3fciifc22)[<S'] 

comparing the second and fourth matrix elements of the left hand side with the second and fourth matrix element of 
the right hand side yields 


^-£^51 {s) 


+ A.is + Bi 

AiS^ -|- A25^ “h A3S “h A4 
4I2S H- B2 

AiS^ -|- A25^ “h A35 “h A4 


Pes2 ( 5 ) 


(A2) 



In the Laplace domain, the turnover time distribution is given by 
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/(s) = k 2 lPESi{s) + k 22 PES 2 {s), 


(A3) 


The inverse Laplace transform of Eq. (IA2I) yield, 


PeSi (t) 


e-‘^\aAi - Bi - a^kii[S]) _ e-^\bAi - Bi - b'^kiijS]) 
(a — b)(a — c)(a — d) (a — b)(b — c)(b — d) 

e-^^(cAi - Bi - c^knlSJ) _ e-^*(dAi - Bi - d^kn[S]) 
(a — c)(c — b)(c — d) (a — d)(d — b)(d — c) 


(A4) 


PeS 2 (i) 


e-‘^*(aA2-B2) e-^\bA2 - B 2 ) 

{a — b){a — c){a — d) {a — b){b — c){b — d) 

e-^\cA2 - B 2 ) _ e-^\dA2 - B 2 ) 

{a — c)(c — 6 )(c — d) (a — d){d — b){d — c) 


(A5) 


where a, 6 , c and d are the effective rate constants, which can be obtained from the numerical solutions of the quartic 
equation + Ais^ + A 2 S^ + A 3 S + A 4 = 0. 


Appendix B: Exact calculation of the turnover time distribution for off-pathway MM mechanism 


The turnover time distribution for the off-pathway MM mechanism [Eq. (1151) ] is given by f{t) = k 2 PESi{t)- This can 
be evaluated from Eq. ([ 6 ]) by taking k[i = k[, k-u = fc_i, fei = ^2 and k [2 = k -12 = ^22 = 0. The resulting coupled 
differential equation, when solved in the Laplace domain, yield the following solution for PeSi ■ 


PeSi (s) 


_-I- B[) _ 

(s^ -I- Ai's^ -I- A 2 's + A 3 ') 


(Bl) 


where, k[ = ki[S], A{ = fci[S'], B{ = ki[S]l3, A/ = fc_i -b fci[S'] -b fe + 2/3, A 2 ' = ki[S]k 2 + A:-i/3 -b 2ki[S]f5 + fe/S and 
A 3 ' = ki[S]k 2 P. In the Laplace domain, the turnover time distribution and its first and second moments are given by 
F{s) = k 2 PESi{s), {t) = —dF{s)/ds\s=o and (f^) = d‘^F{s)/ds^\s=o- The last two expressions yield the mean turnover 
time and randomness parameter [Eig. 5]. In the time domain, on the other hand, the turnover time distribution is given 
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fit) 


[ e-‘^'\-aA[+B[) e-^'*{A[b' - B[) 

^ _ (a' — b'){a' — c') (a' — b'){b' — d) 

e-'^'\-B[+ A'^d) 

(a' — d)(—b' + d) 


(B2) 


where a',b' and d are the effective rate constants, which can be obtained from the numerical solutions of the cubic 
equation + Ai's^ + A2^s + A3' = 0. 
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